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Toward Generalized Entropy Composition 
with Different q Indices and H-Theorem 
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An attempt is made to construct composable composite entropy with diflFerent q indices of 
subsystems and address the H-theorem problem of the composite system. Though the H- 
theorem does not hold in general situations, it is shown that some composite entropies do not 
decrease in time in near-equilibrium states and factorized states with negligibly weak interaction 
between the subsystems. 
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§1. Introduction 

For the recent decade, generalization of the Boltzmann-Gibbs statistical mechanics has gradually 
attracted much attention. Using the (^-modified entropy, so called Tsallis entropyi^ 
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1 



(1) 



many anomalous phenomena, including multifractals, self-gravitating systems, pure electron plasma 
respectively, may be described@'& in "equilibrium" physics language. The index g is a parameter 
unknown a priori and widely believed to be fixed by dynamical details beyond thermodynamical 
feature of the systems. Here the Tsallis probability variables pi are related with the physical 
probability Pi as follows. 

The entropy is also expressed in terms of the physical variables such that 
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Expectation value of energy E is just given in the standard forniB' as 

{E)=Y,E,P,. (4) 



Though recently a subtlety is pointed out when Pj is treated as the fundamental variableJlM* we in 
this paper focus on formulation by use of the variable Pi. It is also possible to discuss a composite 
system of two subsystems A and B with the same q indices. The entropy is just extended like 
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(5) 
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and shows the following non-extensivity for factorized probability distribution Pij = P^P^ , where 
P/^{P^) is probability for the subsystem A{B). 

Sa+b = Sa + Sb + {I - q)SASB- (6) 

Based on the modified entropy one can construct the canonical statistical theory and compare its 
predictions with experimental data of some exotic systems. It is known, for examplejl* that peculiar 
distribution in the quasi-equilibrium state of the electron plasma is well reproduced by the theory 
with q ~ 0.5. 

Although the generalized statistical theory seems to succeeded in phenomenological descriptions 
of many exotic systems in which the standard statistics does not work, we must say that their 
theoretical foundation has not been fully understood. For example, it is still unclear whether the 
second law of the thermodynamics of composite systems with different q subindices holds or not. 
Moreover it is impossible so far to discuss composite entropy of a Tsallis system with q 1 and 
a thermometer made of ordinary matter (q = 1). The essential reason of the poor status comes 
from a fact that guiding principles for determination of the composite entropy have not yet been 
established. 

In this paper, as one of possibilities, we investigate composable composite entropy of two Tsallis 
subsystems A and B with different indices qA and qs- It is shown that bi-composability, especially, 
may possess a powerful ability for restriction of the entropy form. Here the bi-composability means 
that two kinds of composability are simultaneously satisfied. One of the composabilities is just an 
ordinary onefii'i^ 

SA+B{P^3 = P^^Pf) = X{Sa, Sb). (7) 

Thus the total entropy for a factorized distribution Pij = P/^P^ must be computed using only the 
values of subentropies Sa and Sb- By virtue of the composability we are able to take a variation 
of the composite entropy directly using the subsystem variations as 

dX dX 
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Hence it is automatically ensured that the variation of the composite entropy vanishes: 

SSa+b = (9) 

under a deviation as 

6Pij = 5{P,^Pf) = Pf5P,^ + P.Hpf, (10) 

where the distribution Pf-{P^) gives extremum of the subentropy Sa{Sb)- This property will be 
helpful in building concept of thermal equilibrium even in the extended theory. Here the important 
point is that even if subsystems are statistically independent, there is no guarantee that the com- 
posablity should hold. In fact the counter examples can be seen in.El'Ei'0) Another composability 
is introduced when one constructs a grand composite system {A + B) + [A + B)' of two composite 
systems {A + B) and {A + B)' . Let us impose the following property on their entropies. 

S{A+B)+{A+B)' = MS(A+B)-,S(A+By)^ (H) 

where A is an arbitrary function and no need to coincide with the functional form of A. This 
implies that the value of the grand entropy is fixed only by information of the composite entropies 
S(^A+B) ^-iid Si^A+B)' ■ The bi-composability is defined by simultaneous realization of the above two 
composabilities and (0). Recall here that the Tsallis entropy satisfies the bi-composability 
when the q indices of the subsystems are the same. Actually the first composability holds, just seen 
in eqn(^). Moreover a simple computation shows that the second composablity is also satisfied as 
follows. 

S{A+B)+{A+By = S(A+B) + S(A+By + (1 " Q)S(A+B)S{A+By ■ (12) 

Meanwhile it may be instructive to note that some generalized entropies satisfy one of the com- 
posabilities (0) and (|ll| ) but the other not, as seen later. Therefore the bi-composableness is not 
automatically induced, even assuming a single composability (0) or (|ll|) and statistical indepen- 
dence of the subsystems. 

In Appendix A, we give the most general form of the bi-composable composite entropy. In later 
discussions, in order to argue as clearly as possible, we often explain by use of an attractive toy 
model, as an example, called Tsallis-type bi-composable entropy which satisfies a simple pseudo- 
additivity relation: 

S(A+B)+(A+By = S(A+B) + S(^A+By + (1 " Q)S(A+B)S(A+By , (13) 

even when the subindices take different values: qA 7^ qB- The property ( p^ ) is apparently a 
straightforward extension of the result (^) in the same q index case. The parameter Q may be 
regarded as grand q index of the composite systems {A + B) and {A + B)' . It will be often noticed 
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in later sections that a lot of properties which a certain Tsallis-type model has remain still valid 
for more extended bi-composable entropies. 

In this paper we also address the H-theorem problem. Though several works on the H-theorem 
for Tsallis statistics have been already performed,!' §0* the analysis for composite systems with 
different q indices has not yet been discussed anywhere. It is proven in Section 3 that there exist 
some probability distributions in which time-derivative of any composable entropy which satisfies 
eqn (|^) takes negative value. Meanwhile, it is also analytically shown that the H-theorem for any 
state near equilibrium exactly holds for some bi-composable entropies. Moreover, for factorized 
probability distributions with interaction between subsystems negligibly weak, the entropies do not 
decrease in time. 

Noncomposable forms of the composite entropy may be also the other interesting alternative, but 
it looks fairy difficult so far to deal with the entropy systematically, just in a technical sense. It is 
beyond scope of this paper and we expect progress of the investigation in the near future. 

§2. Composability and Composite Entropy 

Let us start to write down the most general entropy form for a composite system A + B oi 
subsystems A and B. The probability is denoted by Pij and the subscript i and j are associated 
with states of A and B. The number of the states of A{B) is denoted by N{M) in the following. 

Note here that for the subsystem A(B) we must respect permutation symmetry in the subscripts 
for the system A{B). Even if the numbering of the A{B) states are rearranged in an arbitrary 
way, the entropy must reproduce the same value as before. Hence the general form reads 



Sa+b — So{Rf, Rf, R2, R2, R3, R3, ■ ■ ■) 



B 



(14) 



where 

M / N \ 

j=i \i=i J 

N ( M \ 

i=i \j=i I 

and 0^, V'^, </)^ and '\\)^ are arbitrary functions. Here is a label of function. 

Next we construct the most general composable entropy which is reduced into a function of two 
subsystem Tsallis entropies Sa and Sb defined as 



(15) 
(16) 



Sa 



Sb 




-QB- 



(17) 
(18) 
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VA = 
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TB = 





Here qAils) denotes the q index of the subsystem A{B). In order to construct the entropy from 
eqn (|l^), it is worthwhile to point out that Tsahis entropies of the subsystems have probabihty 
dependence only through 

Ca = Y.iP^^y' = [1 + (1 - qA)SAV^- , (19) 
i=l 
M 

Cb = Y^iPf y^ = [1 + (1 - <1b)SbV'^ , (20) 

where 

^ (21) 

(22) 

HB 

Taking account of this fact it is easily noticed that the first composability (|^) requires that 

(N M \ 

CA,CB,Y.P^^^llPF\^ (23) 
i=l j=l ) 

where Si is some function and ^ Pf" and ^ will be one after substitution of physical probabil- 
ities because of unitarity. The relation ( p3| ) gives a rather strong restriction to the functions 
and R^. Let us explain the situation using the i?^ case. Substitution of Pij = P-^Pf into R-^ 
yields 

M / N \ I N M \ 

E < E < {p^pf) ]=Ri{cA,CB,Y: pt, T.pn- (24) 

j=l \i=l ) \ i=\ i=l / 

Because there appears a sum over the subscript j on the left-hand-side edge of the left-hand-side 
term in eqn (H), R^ must be always a linear superposition of 

M / N N \ 

Y.ipFy''flAll^p^^y^^llp^^ (25) 

j=l \i=l i=l ) 

and 

M f N N \ 

E^''/21. E(^^r^E^^ • (26) 

j=l \i=l i=l ) 

Here functions f^^ and ^2^1 ^ot determined yet. It should be here stressed again that the partial 
probabilities P^ and P^ must enter into the expressions ( p5|) and (2€) as a simple product P^P^ . 
Therefore the form of R^ is severely constrained and we argue that there exist essentially only the 
following four candidates to take. From eqn (^) the following two components come out: 



M / N N 

E(^''rME(i^^)^- 
j=i \i=i / 
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M / N 

j=i \i=i 



LB. 

"■A 



M / N \ — 

E E(^^.r" -^1' (27) 



j=i \i=i J 

M / N \ '■s 



M / N \ ' 

E(^^r"te^^ 

j=i \i=i / 

M / N \ ' 

E E^/'^" 

7 = 1 \i=l / 



M / N \ '"s 

E E^d -^2. (28) 



j=l \i=l / 

In the same way the following components can be found using eqn (| 
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M / N 

E^^qE(^-'r" 
j=i \i=i / 

M / N \TX 

= E E(^^^^r^ 
j=i \i=i J 

M / N \VX 

= E E(^^.r" ^^3, (29) 
j=i \i=i J 

M / N \ 

j=i \i=i / 

M N 

= EE^^. = i- (30) 
j=i i=i 

The last one is trivial and we have finally three-type building blocks Xa {a = 1,2,3). The same 
argument is also possible for Rj^ and we conclude that the most general form of the composable 
composite entropy is given as a function Q in terms of six independent building blocks Xa and X^ 
as follows. 

SA+B = ^iXa,Xb;rA,rB), (31) 

where 

M / N \ — 

^i = E E(^^.r" ' (32) 
j=i \i=i ) 

^2 = E E^^. ' (33) 

j=l \i=l ) 
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M / N \ — 

^3 = E E(^^.r" ' (34) 
j=i \i=i / 

— 

(35) 
(36) 

(37) 

Here permutation symmetry between the subsystems A and B should be requested in $7, that is, 

0(Xfe, Xa^VB^rA) = n{Xa,Xb; rA, tb). (38) 

So far we have imposed only the first composability ^ on the entropy. It is an interesting 
problem to look for the most general bi-composable form which satisfies the two composabilities 
(^) and (|Tl|). The analysis has been achieved in Appendix A and the result is as follows. 

SA+B = F{A;rA,rB), (39) 

where F{x;rA,rB) is an arbitrary function and 

3 3 

A=f[{Xa)-'"^f[{Xbr''>0, (40) 

a=l fe=l 

and the exponents i^a and v'b are arbitrary constants. Apparently a lot of entropies which are com- 
posable but not bi-composable can be found. Actually it is possible to construct some entropies 
with the form ( pl|) which dependence of Xa and Xi, occurs not only through A. This simple fact 
stresses again that the notion of bi-composability is completely independent of a single compos- 
ability. Moreover if one assumes the pseudo-additivity (|l3|) for the entropy (|39[), it is proven that 
the entropy must take the following form. 

s.,,B=-ie|. (41) 

The proof is given also in Appendix A. The form ( pT] ) is called Tsallis-type bi-composable in this 
paper because double composabilities (0) and (|l|) are simultaneously satisfied and A in eqn (11) 



shows the Tsallis-type non-extensivity with its index Q. The Tsallis-type entropy is just a toy 
model, but as one of the simplest examples it is quite useful to draw what can happen in the 
composing procedure. 

Dynamical details beyond the macroscopic features of the system are expected to determine, in 
principle, the functional form J7, or F, or the values of the indices qa, Qb (and Q for the Tsallis-type 
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model), and the exponents z/^, of the bi-composable entropy, just Uke determination of the index 
q in the case with = qs = q ■ So far we do not know the theoretical method to compute them 
explicitly. However, once one discovers certain composite systems which are described by canonical 
ensembles based on the entropies, it is evidently possible that they are measured by parameter 
fitting analysis of Pij. The canonical ensemble will be given by extremization of the action S 
defined as 



S = Sa+b - a E E ~ 1 




(TV M \ 
EE^^.^^.-(^) ' (42) 
i=l j=l J 

where a and /3 are Lagrange multipliers and generates the unitarity condition: 

N M 

EE^^. = 1' (43) 

i=lj=l 

and the total energy constraint of canonical ensemble: 

N M 

T.T.E^JP^J = {E), (44) 
i=l j=l 

where (E) is expectation value of the total energy and Eij is energy of (ij) state. For example, for 
the Tsallis-type bi-composable case, the equilibrium equation from eqn (E^) just reads 



A 



Q-1 



A Ug dXg ^ 



a + pEij. (45) 



f^XadPij ^^^XbdPij 

The eqns (^) , (^^ and (|45|) generate the probability distribution Pij as a function of qA, qB, Q, 
Vg and i>h- Consequently one can find out the best fit values of the indices and the exponents from 
the experimental data which might discover in the future. 

§3. H-theorem 

In this section let us address the H-theorem problem for the composable entropies. Here we 
assume the positivity of the subindices: ta > 0,rB > 0. Time evolution of the entropy is now 
calculated by using the master equation: 

dP- ^ 

^ = E T.^'^^rMPkl - TM,^JP^j), (46) 
""^ k=l 1=1 

where Tij-^i is transition rate from {kl) state to {ij) state, thus takes positive value by definition. 

Unfortunately the H-theorem does not hold in general situations. It is noticed after rather tedious 
calculation that there exist some probability configurations, governed by a special master equational 
dynamics, in which the whole composable entropy ( |3l[) takes negative values of its time derivative 
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as long as Qa Qb- Thus the H-theorem is actually broken around the configurations. In order to 
explain how it fails, let us write down time derivative of the composable entropy (^), assuming 
the detailed balance situation: 



r 



ij;kl 



kl\iji 



and using the master equation (|4^) as follows. 

dSA+B 



dt 



^ ^ ^ij;kl^ij;kh 
ijkl 



(47) 



(48) 



where 



Sii: 



ij;kl 



E 



-Pki) 

dn (dXa dXa 



dXa \ dPii dP, 



kl 



+ E 



89. / dXi, dXb 



dXh \ dP,j dP, 



kl 



(49) 



Pll 


P12 


Pl3 


PlA 




X + e 





yi3 


yiA 


P21 


P22 


P23 


-P24 







X — e 


y23 


2/24 


P31 


P32 


P33 






y3i 


y32 


y33 


y34 


Pai 


Pa2 


P43 


P44 




Vai 


yA2 


yA3 


I -2x -Y.ycd 



Next let us calculate the 5ii;22 component in eqn ( ^91) for the following configuration in the case 
with (iV,M) = (4,4). 



(50) 



What is firstly noticed here is that 5'ii;22 vanishes if we take exactly e = 0. As the result, expansion 
of 5*11:22 in terms of e behaves as 

Sii,22 = -eK + o{^). (51) 

Consequently if K does not vanish in eqn (|5l|), one can always take negative value of S'ii;22 by 
taking the same sign of the infinitesimal parameter e as that of K: 



5'ii;22 ~ -^K < 0. 



(52) 



Therefore next let us prove the non-vanishingness of K. When we take e = 0, there exist twelve 
independent parameters x and Ucd in the expression of Pij (^0|) . Hence the following twelve variables 
also are functions of the twelve parameters: 



Xa = Xa{x,ycd) (a =1,2,3), 
Xb = Xb{x,ycd) (6 = 1,2,3), 



x,ycd) - T7^{x,ycd) (a = 1,2,3), 



dPii 
dXb 
5Pi 
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cd 



dP22 

dXb 
dp 



22 



{x,y,d) (6 = 1,2,3). 



(53) 
(54) 
(55) 

(56) 
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After a tedious calculation, it can be shown that the Jacobian: 



J 



(57) 



9(x, Vcd) 

does not vanish as long as qA / qs- Thus we can regard Xa-,Xi„Ya and Yf, as free parameters in 
the expression of K. Because the explicit form of K is given as 

^ = E j^{Xa,X,)Y,' + E jr^{Xa,Xtn', (58) 

it is trivial that one can always choose the independent parameters to achieve K ^ Therefore 
we conclude that there exist the configurations for which S'ii;22 is negative. Consequently if we set 
transition rates Tij-ki to zero except 



11;22 — i 22;11 



1, (59) 



decrease of the entropy for the configuration: 

ASa+b 



< (60) 



dt 

is easily shown due to eqn (pSD- Hence we have found breakdown of the H-theorem for any com- 
posable entropy (|3l|) in the configuration (|50|). 

This failure might give a negative impression for the composable entropy itself, but we do not 
think that the situation is so serious. This is because the q- modified statistical picture itself might 
be inadequate to describe the physics for configurations like eqn (^0|). There is a possibility that the 
entropy permits the H-theorem to hold for physically relevant configurations like near-equilibrium 
states, where the q-deformed statistical description is able to work enough. In fact we can find 
some Tsallis-type bi-composable models as the examples in which the H-theorem rigidly holds 
around the equilibrium. When the detailed balance transition rate ( p7| ) is adopted, the equilibrium 
distribution of the model is just equal-weight profile: 

Near the equilibrium the probability is expressed as follows. 

P,j=P,j+eij, (62) 
where e^j are infinitesimal variables and satisfy the unitary condition: 

Eeii = 0. (63) 

For the configuration (^) we are able to compute Sij-ki of the Tsallis-type bi-composable model 
(P) as follows. 

NMA , _ 
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X [{rB{rA - + rA{rB - l)z^i + (r^ - l)i^3 + {^3 - l)i^3}(eii - 



1 ^ 



+{rB{rB - rA)vi + - l)z^2 - (^'yl - 1)^3}-^ ^i^sj - e.sl 



s=l 
M 



+{'rA{rA - rB)vi + r^Cr^ - 1)1^2 - {rB - ^)^3}j^ "^i^js - ^ks) 

s=l 

+o{e^). (64) 

Therefore if one chooses the grand index Q and the exponents fa, z/^ of the model to allow the 
relations: 

rB{rA - 1)1^1 + TAivB - 1)1^1 + (ta - l)z^3 + {rB - l)z^3 „ 
9^ TZTq >0' (65) 

rB^rB - rA)i^i + rB{rB - l)i^2 - {ta - l)z^3 = 0, (66) 
TAirA - rB)i^i + TA^rA - l)z^2 - (?^B - l)i^3 = 0, (67) 

all the coefficients Sij-ki take non-negative values as 

S^r,kl ~ ^gNMA{e^j - e^)' > 0. (68) 



Consequently it is concluded using both eqn (|48|) and positivity of Tij-^i that the value of the 
entropy does not decrease in time around the equilibrium: 

dSA+B ^ Q 
dt 

It is also possible to construct a simple and regular toy model of the Tsallis-type entropy satisfying 
eqns ( |65|) , (|6^) and (^7|) , and it is given in Appendix II. 

Furthermore, beyond the Tsallis-type model, let us introduce a more general bi-composable 
entropy as 

s-^« = o(-ie§) m 

which satisfies 

Ag(x) > 0, (70) 

and eqns (|65|), ( |66|) and (|67|). Then it is trivial that the H-theorem around the equilibrium also 
holds for the entropy (]69| ) due to monotonically increasingness of G. 

It is quite remarkable that the bi-composable entropy ( |69|) does not decrease in time when 
interaction between the subsystems is negligibly weak and initial probability takes factorized form 
as 

P,^{0) = P,^{0)Pf{0). (71) 
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Due to the omission of the interaction, the transition rates can be written as 

^ij;kl = ^i^k^jl + ^ik^f-l, (72) 

where r^^(r^^) is transition rate for the subsystem A{B). The factorized initial condition (71) and 
eqn ( [72[ ) enable the master equation to be reduced into two subsystem equations as 

dP^ ^ A A A A 

= y^(ri:fc-Pfc —^k;iPi )> (73) 
fc=l 

^ = E(r.V."-rf^.^")> (74) 

and guarantees the factorized time-evolution form: 

n,{t) = Pf{t)Pf{t). (75) 

Here it is convenient to introduce a Renyi-type entropy as 

5i? = 3-^lnA. (76) 

Then the time-derivative of the entropy (|69| ) is related with that of the Renyi-type entropy (^) as 
follows. 

dSA^B ^^,(l-^YSn_ ^^^^ 



dt \ I-QJ dt 

Also it is a useful property that the Renyi-type entropy (^) has additivity up to factors like that 

SaiP/^Pf) = caSr^a + cbSr,b, (78) 

where Sji^rA ^^^d Sji^rs the standard Renyi entropies© defined by 



^..^^ln(E(^^)-) 



(79) 



1 



M 



Sr,b = -^ln \j:iPfrn^ (80) 



and CA and cb are constants defined as 



CA = YTT^ ^'^^^^ + ^^^1 + ^^^2 + 1^3) , (81) 
f ~ Qb 

Cb = i_Q K^BVl + TAVl + rBh'2 + i^3) • (82) 

It has been already knowni'i) that the Renyi entropy for positive index monotonically increases by 
the dynamics of the master equation: 

^ > 0, (83) 

^ > 0. (84) 
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By using eqns ( |66|) and ( |67| ) in the expressions of ca and c^, it can be shown that 

CA > 0, (85) 

CB>0 



if eqns (^) holds. Therefore it turns out that time-derivative of the entropy (^9|) entropy takes 
non-negative value: 

dSA+B 



dt 



> 0, (87) 



dt dt ~ 



AG' ( --^ — ^ ) > (89) 



because the inequality: 



holds and 



l-Q 

is guaranteed due to eqn ([70|). Consequently it has been proven that there exists the composable 
entropies which satisfy the H-theorem for both the near-equilibrium case and the uncorrelated 
subsystems case, and an explicite example is the form (|69|). 

It may be interesting to comment, along Abe's argumentlll) (also see , on the thermal balance 
relation generated by the entropy (|69[), though the discussion is still in a rather heuristic level. It 
turns out that the nonextensivity of the entropy (^) is given as follows. 

Sa^b = G (-^) > 

A = [1 + (1 - q^)SAYBVl+v,+rAV,+rAV2 

X [1 + (1 - g^)5^]^^APl+P3+rs-i+r-B-2 _ (90) 



Using the H-theorem conditions (|66| ) and (pTI), the variation of the entropy can be written as 

.c- ^B{rB - l){yi+T^2)+rA{rA - l)(i^i + 2^2) . ^, 
5Sa^b = 

9^ SSa + -—^^——5Sb 



l + {l-qA)SA l + {l-qB)SB 

Then we take 5Sa+b = under the total energy conservation relation: 



(91) 



5Ea + 5Eb = (92) 

to get the thermal balance relation. The result is given as follows, independent of the functional 
form G{x) and the value of Q. 

qA SSa qB SSb 



1 + (1 - qA)SA SEa 1 + (1 - qB)SB SEb ' 



(93) 
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Note that this relation includes the Abe's balance relationEl^ as a special case. Actually when 
QA = Qb = Q taken 

^ ^ = I ^ m) 

l + {l-q)SASEA l + {l-q)SB5EB ^ ^ 

is exactly reproduced. Also eqn ( |93| ) is consistent with a guessed relationll' by Tsallis for the different 



q case. Though we have a subtlety that the canonical thermal state of the composite noninteracting 
system is not the product of the states of subsystems for the generalized entropy, there is a possiblity 
that the above result ( p3D remains unchanged for thermodynamic-limit situations just as in the same 
q caselli) (also seei^ ) . 

If the system A is taken as an ordinary Boltzmann-Gibbs system {q^ = 1), the subentropy is 
reduced into the BG form Sbg-.a- Then the relation (|9^ is re-expressed as 

SSbg-.a Ob SSb , . 

5Ea " 1 + (1 - qB)SB 5Eb ^ ' 

Even in the extotic situation, the fact can be considered robust that physical temperature Tpj^yg 

can be introduced for the system A as follows. 

1 SSbg:A 



Tphys SEa 

Therefore observable temperature Tg of the Tsallis system B should be defined as 

1 _ Qb SSb 

Tb ~ l + {l-qB)SB6EB' 

so as to preserve the zero-th law of thermodynaics: 

1 1 



(96) 



(97) 



(98) 



'Tphys 'Eb 

Here we want to mention that before our analysis no one argues explicitly presence of the numerator 
of the prefactor in the right-hand-side term of eqn (|97[). Due to the definition (|97[), the original 
relation (pSj), in which qA is not needed to take unit, can be interpreted as a generalized thermal 
balance as follows. 

^ = ^- (99) 
-tA J-B 

The transitivity relation (^) looks quite plausible and attractive, though the derivation seems still 
nonunique. 



In summary, though the H-theorem of all the composable entropies ( |3l|) does not hold for some 
probability distributions like in eqn (|50|), it has been rigorously proven that there exists a rather 
general entropy (^9|) which does not decrease in time for both any states around the equipartion 
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equilibrium and the factorized distributions without interaction between the subsystems. We also 
derive the thermal balance relation for the different-indices case from a plausible but still heuristic 
argument. 

Though analysis for noncomposable composite entropy has not been included in this paper, 
the possibility remains, of course, still alive and may be interesting. However the systematic 
investigation seems complicated and seems to need some innovation in the future. 
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Appendix A: Bi-Composable Entropy 

In this appendix, let us prove first that the most general bi-composable entropy takes the following 
form. 

Sa+b = F{A;rA,rB), 

where the function F is arbitrary and A is defined in eqn (|40|). Recall here that the most general 



composable entropy is given in eqn (31). Thus our remaining task is just imposition of another 
composability (11) on J7 in eqn (^T]). To achieve this, it is very helpful to note a fact that 



N N' 

EE 

1=1 i'=i 

N 

E 

,i=i 



(A+B) p{A+By 
ij i'j' 



M M' 

EE(^. 

j=i f=i 

-M f N' 



i:(^? 

i=l 



E 

\i'=l 



M' 



■By 



i'=i 



(A-l) 



Thus all the six components Xa and Xh are completely factorized for the probability forms as 

(A-2) 



p _ p{A+B) p{A+By 
^ijkl — ^ij J^^kl 



Consequently we obtain the following parametric relation among S(^a+b)-, S(^a+b)' ^-nd 



(A+B)+(A+By ■ 



5, 



{A+B)+{A- 



-By = ^{Sata,Sbtb;rA,rB), 



where 



S{A+B) = ^{Sa,Sb;rA,rB); 
S(A+By = 0,{ta,tb;rA,rB) 



(A-3) 
(A-4) 
(A-5) 

(A-6) 
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''a ^'-a 1 

(A+B) 
b ' 



Sb = X, 



(A-7) 
(A-8) 
(A-9) 



can be regarded as free parameters of the expression, independent from each other by taking 



N,N',M and M' largely enough. Using the relations ( |A-3| )'^( A-5 ), it is possible to write down a 
lot of parametric expressions of 5'(A+_B)+(yi+B)' terms of S(^a^b) and 5'(a+b)'- For example, we 
obtain an expression by taking 



Sl = S2 = S3 = Si = S2 = S3 
^1 = ^2 = ^3 = ^1 = ^2 = ^3 = 



as follows. 



S(^A+B)+{A+By = ^{(t>G, (pO, (j)0, (j)0, (pe, (pO; rA, vb), 

S{A+B) = ^i4>, 4>, 4>, 4>, 4>, 4>; rA^rs), 

S^A+BY =^{0,e, 9,6,0, e;rA,rB). 



(A-10) 
(A-ll) 

(A-12) 
(A-13) 
(A-14) 



The above relations ( A-12| )~( A-14| ) indicate that S(^a+b)+{A+by is given as a function of S(^a+b) 
and S^A+B)' implicitly via the two parameters (j) a-nd 0. 

Next we must find explicitly the general solution of eqns (|A^)^([A^. In order to do that, let 
us consider first the following relation obtained from eqns ( [A -41 ) and ( [A-13 ). 



i^isa,si,S2,S3;rA,rB) 

= ^i4>A,4>,4>,4>A;rA,rB)- 



By solving eqn ( A-15 ) for S3 , we can define a function a as 

S3 = Cr{4',Sa,Si,S2). 

Substituting eqn ( A-16| ) into eqn ( A-3 ) and using eqn ( |A-12D yield 

^{sah, siti, S2t2, cr{(t>, Sa, Si,S2)cr{9, ta, hM)'-, rA,rB) 
= n{ct>0,(P0,(P0,(t>0,(t>0,(t>0;rA,rB). 

Then it is soon noticed that eqn ( A-17| ) means 

(j{4), Sa, Si,S2)cr{9, ta, hM) 
= Z{Sata,Siti,S2i2,' 



(A-15) 



(A.16) 



(A-17) 



(A.18) 



where Z is & certain function. Eqn (|A-1^ ) looks very complicated, but fortunately can be solved 
by treating the variables as six pairs: (s^, ia),(si, ti), (s2,t2) and {4>,9). For each pair we solve the 
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equation independently. Actually for an arbitrary pair {ps,pt) in the six, eqn ( A-18 ) is regarded 
as 



(A-19) 



^From eqn ( A-19| ) the following relations are straightforwardly obtained by taking ps = 1 or pt = 1. 

1 



Hp) 



a(i) 
1 



Zip), 
Z{p). 



Substituting eqns ( [A^ and (lA^) , eqn ( [A- 191) is replaced into 

Z{ps)Z{pt) = Z{l)Z{p,pt). 



(A-20) 
(A-21) 

(A-22) 



By taking differentiation with respect to pt in eqn ( A-22 ) and setting pt = 1 and ps = p, the 
following equation appears. 

dZ _ Z'{1) 



p-r 



-Z{p), 



(A-23) 



dp Z{1) 

where Z'{1) is derivative of Z at p = 1. This equation is easily integrated and leads us to the 
following relation: 

1 



HI 



■Z{p) OC /3S 



(A-24) 



where ^ is a constant. After the same procedure has been performed for all the six pairs, it is finally 
found that 



S3 = <y{<P,Sa-,Si,S2) OC 



(A-25) 

a=l 6=1 

where the exponents £,^,S,a and ^j, are constants. Solving eqn ( [A -251 ) for (j) leads to the following 
result. 

3 3 

ct> = Co\{{sa)-''-\{{sb)-'' = Co/^, (A-26) 

a=l 6=1 

where Co, i^a and Vh are also constants. By introducing a new function F as 

F{x;rA,rB) = 0,{CoX,CoX,CoX,CoX,CoX,CoX;rA,rB), (A-27) 

the final result ( |39|) is achieved using eqns ( [A45D and ( |A^ . 

Next let us prove that the Tsallis-type bi-composable entropy takes the form (|4l]). Firstly let 
us impose the grand pseudo additivity ( [T3D on the bi-composable entropy (p9|). Then it turns out 
that the function F possesses the following property: 



F{xy; rA, tb) = F{x; ta, tb) + F{y; ta^vb) 

+ (1 - Q)F{x; rA, rB)F{y; r^, vb) 



(A-28) 
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for arbitrary parameters x and y. Differentiating eqn ( A-28| ) with respect to y and taking y = 1 
yield 

dF 



x-^ = F'{1; ta^tb) + (1 - Q)F'{1; rA,rB)F{x; r^, r^), 



(A-29) 



where F'{l;rA,f'B) is derivative of at x = 1. Eqn ( A-29 ) can be integrated and we get the 
following solution: 

1 



F{x;rA,rB) 



l-Q 



+ C'x'^ 



(A-30) 



where C and rj are constants. Substituting eqn ( A-30| ) into eqn ( [A-28 ) fix the constant C as 

1 



c 



and finally the function F is determined as 

F{x]rA,rB) 

Therefore, redefining as 



l-Q 



l-Q 



[1 - x''] . 



(A-31) 

(A-32) 

(A-33) 
(A-34) 



in the definition of A, it is concluded that the most general form of the Tsallis-type bi-composable 
entropy is given by eqn (41). 

Appendix B: A Regular Example of Tsallis-Type Bi-Composable Toy Model 

Let us discuss a plain case that the grand index Q and the exponents Ua, Vb behave as func- 
tions only in terms of va and tb- Then in the functions the permutation symmetry between the 
subsystems A and B should be preserved: 



Q{rA,rB) = Q{rB,rA), 

h{rA,rB) = i^b{rB,rA), (6=1,2,3). 



(B-1) 
(B-2) 



Furthermore when we take = tb, that is, qa = Qb, the composite entropy must be reduced into 
the original Tsallis entropy (|5|). Hence we require the following boundary conditions: 

1 



Q(r,r) = - = q, 
r 

I . I q 

^ ' ' 2r 2' 
2^2 (r,r) = z^3(r, r) = 0. 

Then let us introduce a entorpy form as 

rA{rA-l)^ + rB{rB - if 



(B-3) 

(B-4) 
(B-5) 



Sa+b 



{rA-l){rB -l)irA + rB -2) 
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1 - xl^-'xl 



1 













(B-6) 



Clearly if = fB is taken, the entropy in eqn ( |B-6| ) is reduced into the original Tsallis entropy. 



Also it can be pointed out that eqns (pq), (pq) and (pTI ) certainly hold in eqn ( B-6 ), thus the entropy 
satisfies the H-theorem for near equilibrium and factorized states. 
Here one might worry about a "singularity" at 



rA + rB = 'i 



(B-7) 



in eqn ( [B-6D . But fortunately it is explicitly checked that if one takes — > 2 — limit, any 
singularity does not take place and the entropy really has well-defined finite form as 

2\ 



hm Sa-\-b 



2(1 -ta) \Xi \X2) I' 
The limit of r^i — > 1 is also obtained as a regular form. 



where 



X 



l;2 



lim Sa+b 



Xo T H : r 

2X2 2 ^'^ 2(1 -rs) 



^ In X2 + In X3, 



2(1 -tb) 



lim 



Xi — X2 



rA^i — 1 
M / N 



M / N \'^B /AT \ M / N \ ^B~l / TV 

E^d In E^^'. +-bE E^d E 

7 = 1 \j=l / V = l / 7 = 1 \j=l / \i' = l 



X2-I = lim 



X2-I 



r-A^l rA — I 
N / M \ / M 

E E^din E^^.' 



fB- 



(B-9) 
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